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Abstract 

We analyze the scalar radiation emitted from a source rotating around a 
Schwarzschild black hole using the framework of quantum field theory at the 
tree level. We show that for relativistic circular orbits the emitted power 
is about 20% to 30% smaller than what would be obtained in Minkowski 
spacetime. We also show that most of the emitted energy escapes to infinity. 
Our formalism can readily be adapted to investigate similar processes. 
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I. INTRODUCTION 



Observational confirmation of the existence of black holes is one of the most important 
challenges in astrophysics. Recently a number of compact objects in X-ray binary systems 
have been identified as black holes since a careful analysis has shown that their masses are far 
beyond any limit accepted for dead stars in general relativity Jl|. There also exists indirect 
evidence of the presence of supermassive black holes in the centre of some galaxies 0. 
Nevertheless, unambiguous confirmation of the existence of black holes would require the 
observation of effects due to the event horizon itself. This is expected to be achieved through 
precise measurements of the electromagnetic radiation emitted from black hole accretion 
disks (see e.g. and from the gravitational radiation emitted from companion stars 

orbiting black holes (see e.g. 0). Because radiation from sources orbiting black holes plays 
such a crucial role in modern astrophysics and because increasingly precise measurements 
are leading to the observation of relativistic effects occurring in the vicinity of the horizon || , 
investigation of how radiation-emission processes are modified by the nontrivial curvature 
and topology of the black- hole spacetime is particularly important. 

In this paper we analyze analytically and numerically the scalar radiation emitted by a 
source rotating around a black hole using the framework of quantum field theory at the tree 
level, and compare the results with those obtained in Newtonian gravity and in a theory 
associated with one-graviton exchange in flat spacetime. We show that these results coincide 
asymptotically, as expected, but considerably differ close to the last stable circular orbit. 
We also calculate the amount of emitted energy which is not absorbed by the black hole. 
The paper is organized as follows. In section || we present the framework of quantum field 
theory in which we will work. In section [111] we obtain analytical and numerical results 
for the emitted power. In section |y| we establish a connection between our quantum field 
theory approach and classical field theory. In section |V| we compare our curved spacetime 
calculation with flat spacetime calculations by considering (i) Newtonian gravity and (ii) 
the one-graviton exchange theory. In Section |Vj we calculate analytically and numerically 
the amount of radiation which is not absorbed by the hole. Finally in section |V11| we make 

some remarks on our results. We use natural units c — h — G—1 and signature (H ) 

throughout this paper. 



II. GENERAL FRAMEWORK 

A nonrotating black hole with mass M is described by the Schwarzschild line element 

ds 2 = f(r)dt 2 - fir-Y^r 2 - r 2 d6 2 - r 2 sin 2 9 d(f) 2 , (1) 

where f(r) = 1 — 2M/r. Let us consider a circularly moving scalar source at r = Rs with 
constant angular velocity Q > (as measured by asymptotic static observers) on the plane 
9 = 7r/2 described by 

j'M = 7=v^ r ~ Rs) 6{e ~~ n/2) 5(0 ~ nt) (2) 

with g = det(g^ u ), where the constant q determines the magnitude of the source-field cou- 
pling. The four-velocity of this source is 
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Rs) = [(f(R s ) - Rln 2 )- 1 / 2 , o, o, n/(f(R s ) - R 2 n 2 y 



/21 



(3) 



We have normalized the source j by requiring that / daj{x^) = q, where da is the proper 
three- volume element orthogonal to u^. Let us now minimally couple j{x^) to a massless 
scalar field ^(x^) so that the total Lagrangian density is 



(4) 



The positive-frequency modes in spherically symmetric static spacetime can be given in 
the form 



7r r 



Yr 



Irri \ 



-iuit 



{to > 0) 



(5) 



where Ou^im = 0. In the present case the functions V^zfV) satisfy the differential equation 

VS(r)=wfyS(r) (6) 
with the following scattering potential (see figure IT]): 



-/W^f/(r)^)+^(r) 
ar V dr I 



V s (r) = (1 - 2M/r) 2M/r 3 + /(/ + l)/r 



(7) 



In terms of the dimensionless tortoise coordinate x = r/2M + ln(r/2M — 1), can be 
rewritten as 



d 2 



dx 2 



+ AM 2 V s [r{x) 



(8) 



The scalar field can be expanded in terms of the complete set of positive- and negative- 
frequency modes as 



$ in (^ 



oo I 

E E 



du 



o 



1=0 m=-l ' 

The Klein-Gordon inner product a kg is defined by 



(9) 



a KG {<P, V) = % I dL n^(0* - V^0* • V) 
7e 4 



(10) 



where n^ 1 is the future-pointing unit vector orthogonal to the Cauchy surface S t with £ 
const. The modes « w ; m are Klein-Gordon orthonormalized: 



VKGiUulrm U^'l'm') — S(uJ - Uj')Su>S ri 



(11) 

(12) 



Then the creation and annihilation operators, a^i m and a™ lw satisfy the usual commutation 
relations [a™ lm , dj Vm ] = 8{uj - u')6u>6 mm r. 
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We see from (RI) that close to and far away from the horizon, the modes ip^ (x) purely 

incoming from the past horizon H~ , and the modes ip^i (x) purely incoming from the past 
null infinity J~ can be written as [|7].|j 

Xs (x) w I A ujl (e 2 ^+ jLi e~ 2tM " x ) (r > 2M) , 
u \ 2i l+l A wl Tui Mux h\ 1] (2Mujx) (r > 2M) , 

and 

A (x) * I Kl (r - 2M) ' _ (14) 

| 2A ujl Mujx[(-i) l+1 hf ) (2Mujx) + TZ^i (2Mlux)} (r»2M) , 

where h\ (i = 1, 2) are spherical Bessel functions of the third kind || and the overall 
normalization constant is determined by ( |TT]) as = (2a;) -1 . 
The emitted power is given by 

W t Z= l +0 ° duu\A^ m \ 2 /T , (15) 
Jo 

where 

= (ulm\ i I d A x^9 |o> = < / dW^i^R/J^) ( 16 ) 



is the emission amplitude at the tree level, and T = 2n5(0) is the total time measured by 
asymptotic observers |10| . We have chosen the initial state |0) to be the Boulware vacuum, 
i.e. a^ m |0) = 0. If we had chosen the Unruh or Hartle-Hawking vacuum flllj, then (15) 



would be associated with the net radiation emitted from the source since the absorption 
and stimulated emission rates (which are induced by the presence of thermal fluxes) are 
exactly the same. Note that for sources in constant circular motion the amplitude A e ^ l rn is 
proportional to 8{uj — mVl). Hence the frequency of emitted waves is constrained by us = mQ. 
In particular, since Q > 0, no waves with m < are emitted. 



III. EMITTED POWER: NUMERICAL AND ANALYTICAL RESULTS 



The general solution of (HD is not easy to analyze [12] and therefore we first esti- 
mate the radiated power numerically. (Later we will compare it with the one obtained 
with our analytic approach.) To do so we solve @ for the left- and right-moving ra- 
dial functions ip^i ( r ) an d ip^i ( r ) with asymptotic boundary conditions compatible with 
equations (|IB|) and (|T4T). The method for finding the the functions ip^i can be sum- 
marized as follows. (The functions ip^ are obtained in a similar manner.) We recall 

from (|14D that close to the horizon ip^cc exp(— 2iMux). Thus we construct the solu- 
tion Xui( x ) °f ® with fixed to satisfying Xui{ x ) = exp(— 2iMux) as the initial condi- 
tion for x < 0,\x\ ^> 1 by evolving it numerically towards large x. Far away from 
the horizon we have Xui{ x ) ~ Aexp(— 2iMujx) + B exp(2iMux) [see ( |i~4"|) and recall 
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that hf ) (2Mujx) ss (-i) l+l e 2lMuix / (2Mujx) and hf ) (2Mujx) ss i I+1 e- 2iAfl "7(2Mu;x) for 
2Mcja; ^> 1]. The constants A and 5 (with |v4| 2 — \B\ 2 = 1) are determined from the 
numerically obtained values of Xwi( x ) and dx^i/dx for 2Mux ^> 1. Now, by requiring the 

asymptotic boundary conditions compatible with (|T3J), we find ^ (x) = Xwz(^)/(2o;A). (In 
particular, T^= 1/M. and 7^= B /A.) 

By substituting the radial functions t/q and ^ in (g) to construct the right- and left- 
moving modes u w j m and u u i m , respectively, and using (|j"5|) with (|16[), we find the corre- 
sponding radiated powers: 

WC= 2m 2 n 2 q 2 [f(R s ) - R 2 s tf} | ^ (i? 5 )ARs| 2 1^/2, m)| 2 , (17) 

^S,em 

where uq = mO and /,m > 1, and similarly for W\ m • We note [IB] that Yi m (7r /2,Qt) = 
if I + m is odd and 

W UlOOtV* 2l + l (l + m-l)\\(l-m-l)\\ 

lYU^m (l + m)\\(l-m)\\ (18) 

if I + m is even, which is of course time independent. We have defined n\\ = n(n — 2) ■ • • 1 
if n is odd and n\\ = n(n — 2) • • • 2 if n is even and (—1)!! = 1. Now, the condition that our 
source be in a stable circular geodesic implies Rs = (MVL~ 2 ) 1 ^ as is well known ||14|| . By 



> em <— em 



using this formula, Wi m and Wi m can be cast as functions of quantities measured at infinity 
alone, namely, Q and M. In figure |2] we plot the total radiated power with fixed angular 

S em ^S,em <—S,em 

momentum, W lr ^ m = Wi m + Wi m , as a function of the angular velocity f2 for different 
values of I and m. Note that W^ m = for odd / + m because Yj m (7r/2, Vtt) = in this case. 

Next we consider an analytic approximation valid for low-frequency modes. Let us first 
recall that waves emitted from circularly moving sources obey the constraint u = mQ. 
Thus, waves with m = 1 (which turn out to be the most important ones for our purposes) 



emitted from sources in stable circular geodesic orbits (Rs > QM, Q = JM/R S ) have 



maximum frequency aj max = (Qy/Q M) 1 <C y Vg" ax , where is the maximum of the 

scattering potential Vs. Hence waves with small angular momentum have small frequencies 
in comparison with yfv^, i.e., w 2 /V 5 max < (u max ) 2 /Vg 13 * « 4 x 1(T 2 < 1 (see figure @). A 
similar analysis for waves with arbitrary m shows that u 2 /V™* < 10 _1 . Therefore, we will 
consider the approximation where the radial functions are replaced by their leading terms 
for small uo (see equations (6.6) and (7.2) in ||, and |L5| for a correction): 



^ l (r)^2rQ l (r/M-l) (19) 

and 

Ts (r) „ 2 2l (l\ f(Mu J yrP l (r/M-l) 

^ {r} (20! (21 + 1)! ' 1 Uj 

where Pi(x) and Qi(x) are the Legendre functions. By substituting ip^ and ^ (given in 

(-em 

equations (|19"D and (p0|) respectively) in ([[7]) and a similar expression for W\ m , we find the 
corresponding radiated powers [Rs = (Mfi -2 ) 1 ' 3 ]: 
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8q 2 m 2 Q 2 (f(R s ) - R 2 S Q 2 ) \Qi(R S /M - 1)| 2 \Y lm (n/2, Qt)\ 2 (21) 



and 

(f(R s )-R 2 s n 2 )\Pi(R s /M-l)\ 2 \Y lm (*/2,m 2 (22) 



-s,em^ 2 il+1 q 2 (l\) 6 m 2l+2 M 2l n 2l+2 
Tlm ~ [(2/)!] 2 [(2/ + l)!] 2 



where l,m> 1. The total radiated power with fixed angular momentum obtained in this 
approximation is also shown in figure ^. We see from this figure that our analytic approx- 
imation has better accuracy for small Q as expected. [Note that u oc Q and see discussion 
above (0)]. 

It can be seen from figure |2| that the total radiated power 

oo I 

w S,em = £ £ ^ (23) 
Z=l m=l 

will be dominated by waves with small I. For circular geodesic orbits far enough from the 
horizon we can use equations (0)-(23) to write the radiated power in a simple form: 



W 



S.em 



R>>rs ^q 2 M 2 ^n 8 / 3 /12n . (24) 



IV. CONNECTION WITH CLASSICAL FIELD THEORY 

Since our calculations are performed at the tree level, our results can be interpreted in 
classical field theory as follows (see, e.g., fl6|). One can show that the energy of a classical 
field can be written as 

E= l -a KG ^,d t <P) , (25) 

where the Klein-Gordon inner product a kg is defined by fllPl) . The classical field generated 
by a source j{x^) can be expressed in general as 



(a?) = % 



EE / du [A e j m u.uxn - Atz^uxn*) , (26) 

;=o m=-l J0 

where A e l ^ m is defined by equation (|T6|) . By substituting this in (p5|) and using the orthonor- 
mality of u^im we obtain the energy as 

E = E E / u;\A e J m \ 2 , (27) 



1=0 m=-l 







which agrees with fll5|). (See, e.g., |17j for other classical- field-theory analyses mainly devel- 



oped to study gravitational wave emission.) 
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V. CURVED VS. FLAT SPACETIME CALCULATIONS 



Now we compare the radiated power W ,em in Schwarzschild spacetime with that in 
Minkowski spacetime. We consider radiation from a source in circular motion in Minkowski 
spacetime due to the presence of some gravitational force. The latter should give fairly good 
results for the case of a source rotating around a star that is not very dense, but not for the 
case of a source rotating close to a black hole as will be shown. 

Let us represent the scalar source in Minkowski spacetime by 

/V) = -J— S(r - R M ) 6(6 - tt/2) 6(4> - Qt) , (28) 

K M 7 

where 7 = 1/ ^Jl — R 2 M VL 2 . We are working here with spherical coordinates defined through 
the Minkowski line element given by ([[]) with f(r) = 1. Note that the normalization of 
this source is chosen so that / da j M (x^) = q, where da is the proper three- volume element 
orthogonal to the world line of the source, as in the Schwarzschild case. 

It is not a priori meaningful to compare the radiated powers from the sources in 
Schwarzschild and Minkowski spacetimes with the same value of r. For this reason, our 
results will be expressed in terms of Q. It is meaningful to compare the radiated pow- 
ers using Q because this is a (coordinate free) quantity measured by asymptotic observers. 
We expand the scalar field as in equation @ with Klein-Gordon orthonormalized positive- 
frequency modes given by @, where the radial functions are given by V'*[( r ) = r ji( ujr ) 
(see, e.g., [|J), which satisfy the differential equation 

(-^ 2 + M^W=^W (29) 

with Vm = 1(1 + l)/r 2 . The total radiated power calculated in Minkowski spacetime is 

00 1 

W M,em = £ £ 2q 2 m 2 n 2 1 - 2 |jKmfii? M )| 2 |^ m (7r/2,^)| 2 , (30) 

1=1 m=l 

where we have used the Minkowski vacuum as the initial state. In order to compare this 
expression with the radiated power obtained through classical field theory, we adapt the 
standard derivation of Larmor's formula for electric charges (see, e.g., |18[) to the case of 
scalar sources: 

W^T = <A7l2vr , (31) 

where a is the proper acceleration of the source with an arbitrary trajectory. For circular 
orbits, we have = 7 2 n 2 R M . The equality W M ' em = W^T follows from the following 
formula: 

00 1 1 2 

E E m%( m z)f\Y lm (^2A)\ 2 = ^TT^-m ' ( 32 ) 

1=1 m=l Zq:7T \ l Z I 

for \z\ < 1. A proof of this formula is given in the Appendix. We also verified it numerically 
for a wide range of values of z. 
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Now, in order to compare W M,em with W s ' em , we shall cast Rm as a function of Q by im- 
posing the condition that the scalar source be in circular orbit due to the influence of a grav- 
itational force. There is no unique way to define a gravitational field in flat spacetime (see, 
e.g., ||19|| ). Assuming Newtonian gravity and using Kepler's third law Rm(Q) = (Mfi -2 ) 1 / 3 
in (gO]), we find W M ' em as a function of fi. We plot the ratio W s ' em /W M ' em in figure | as a 
function of Q. A similar result is obtained if one considers the gravitational force generated 
by the one-graviton exchange diagram [[TJJ. In this case a straightforward calculation gives 
7 2//3 i?Af = (Mfi" 2 ) 1 / 3 where 7 is defined below equation (|28|). By solving this equation for 
Rm, we obtain 

Rm{M) — — r + 

3 (27 M fi 4 - 2 M 3 fi 6 + J -A M 6 fi 12 + (27 M fi 4 - 2 M 3 fi 6 ) 2 ^ ' 



27 M fi 4 - 2 M 3 fi 6 + J-4 M 6 ft 12 + (27 M fi 4 - 2 M 3 ft 6 ) 



M 



+ 1 , (33) 

2^3 ft 2 3 1 ; 

where R M > 0. We substitute (|3]) in (|30|) to find the radiated power W M ' ern . We plot in 
figure [| the ratio W ,em /W M,em as a function of ft. Asymptotically, the radiated power is 



W 



M,em 



2 M 2/3 fi 8/3 /127r (34) 

in either case. Figure [3] and figure ^ show that for ft — > the ratio W s,em /W M,em tends to 



unity as it should [see equations fl2~4f ) and @)]. As the source approaches the last stable orbit, 
however, our numerical and approximate analytic results show that W s,em /W M,em < 1. This 
is not a trivial consequence of redshift since the frequency ft is measured at infinity in both 
the Schwarzschild and Minkowski calculations. 



VI. ASYMPTOTIC RADIATION 



Finally, it is interesting to compute the amount of the radiated power that escapes to 
infinity, namely 



00 1 

w s,obs = y^J2 

1=1 m=l 



— » ^ ^S,em <— ^ *—S,em 
TuqI I Wim +| Kw l I W irn 



(35) 



By using | TZ^i | 2 = | 1Zuj i 1 2 — 1 — I Tu i \ 2 , the observed power (R3) can be written as 



W 



S.obs 



00 I 

EE 

1=1 m=l 



. S.em 



r^ i nw lm - w lm )+ w 



- S,em' 
Im 



(36) 



A numerical estimate of W ' obs /W ,em is given by the solid line in figure [5]. In order to obtain 

an analytic approximation to the transmission coefficient | Tuj i | 2 and hence to W s,obs (for 
uJq/Vs 13 * <C 1) we first note that asymptotically fll9|) reads 



S 



<i (0 - fOT r » 2M • ( 37 ) 



Now, by fitting ( B7| ) with ([13]) (with = cjo) in the range of r satisfying r ^> 2M and 
u r C 1, we obtain 

2 2 '+ 2 (Z!) 3 (M^ o y+ 1 
' °' ' (2/ + l)!(2/)! ' 1 j 

The analytic approximation of W s,obs /W s,em obtained in this way is given by the dashed 
line in figure [5[ It is seen that very little of the emitted radiation is absorbed by the black 
hole. This does not contradict the fact that Schwarzschild black holes have a non-negligible 
absorption cross section for infrared particles p0] (actually of the order of the horizon area), 
because the main contribution to the cross section comes from modes with I = 0, which are 
not emitted by our circularly moving source. 



VII. FINAL REMARKS 



In summary, we have calculated the radiated power from a scalar source rotating around 
a black hole in the framework of quantum field theory at the tree level. We have shown 
that for relativistic circular orbits the emitted power is about 30% and 20% smaller than 
what would be obtained in Newtonian gravity and in flat-spacetime theory with gravitation 
generated by one-graviton exchange, respectively. This is in agreement with the fact that 
astrophysical processes involving wavelengths of the order of the Schwarzschild radius need 
to be described using fully curved spacetime. We have also shown that most of the emitted 
energy escapes to infinity. Clearly, the presence of surrounding matter could trap the emitted 
radiation in the vicinity of the black hole because of friction effects [[§]. These astrophysical 
issues, however, are beyond the scope of the present paper. Our conclusions are qualitatively 
in agreement with those obtained for gravitational waves. The procedure used here can be 
readily adapted for scalar sources following other trajectories (see, e.g., [^TJ for chaotic ones) 
by altering the scalar source (Q). 
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APPENDIX: A PROOF OF EQ. (32) 



In this Appendix we prove the following formula: 

oo I 1 Z 2 

24vr (1 - z 2 ) 3 



/W=EE^bK^)] 2 |^/2,0)| 2 = -^-^ (M<1). (32) 



1=1 m=l 
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We note here that \Yi m (n/2, 0)| 2 = \Yi m (n/2,0)\ 2 is independent of 4>. 
We start by showing that 

/(*) = *T2 ^ f ^ sin b(* sin ^)] 2 ' ( 39 ) 
62tt a Jo Jo 

where 

/+oo 
dip 5' (<f) - ip + a sin ip) , (40) 
-00 

if — 1 < a < 1. By combining the formulae 



oo 

j,kz cos 7 



£(2Z + l)i'j,(A;z)Pi(cos7) (41) 

1=0 



and 



47T ' 

P*(cos T ) = — — 2 Yim(e',i>)*Y lm (8A) , (42) 

m=—l 

where cos 7 = cos 9 cos 9' + sin$sin#'cos(0 — ?/>), and letting 9' = 7r/2, we have 

00 i 

e i fc ,sin co S (^) =47r ^ £ ^(^)F /m (7r/2,^)*F /m (^0) , (43) 

£=0 m=-Z 

By using the formula (27r) _1 g^/> e lM ^ = 5(p), we find 

1 /" + 0O 

/ j^^p ^ikzsinB cos((p—ip)+ikip 



2tt 



00 

00 i 



= 47r E E ^iK^)^m(7r/2,O)^,0)5(A;-m) . (44) 

l=0 m=-l 

By multiplying by k and integrating over we have 

i r+00 r+00 

G(z,9,(f)) = — dk d ^ ke ikzsin8co S (<P-^ +l ki, 



57T 2 J-00 



00 I 

= EE* ,m ^W^m(7T/2,O)[^ ro (0,0)-^ i _ m (0,0)] . (45) 
i=l m=l 

(We are using the convention Y/ i __ m (6',0) = Yj m (0, 0) here.) Notice that this function is 
periodic in with priod 2it. By using orthonormality of the spherical harmonics Yi m (9, <fi) 
we find 

f(z) = - / #/ ^sin#|GM,0)| 2 . (46) 
2 Jo Jo 

By shifting the integration variable in fl45|) as ip — > ip + <f), we have 
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G(z,e, 



dk 



57r J-oo 



d 



dip k e ik ^ z sin e cos ^+^ +< W 



-oo 
oo 



571"° J-oo 
l-OO 



+00 



dip / dke ik{z sin9 ms ^+<t>) 



1 r+°° 

— — / dip 8 (zsine cosip + ip + 

4:71 J-OO 



(47) 



By changing the variable as ip — > tt/2 — ip we find 

1 r+°° 



1 f+°° 

G(z,9,(p) = -— / dip 5'(<f> + tt/2 - ip + z sine sin ip) 

4lT J-oo 



-~T9{z sine, + 71-/2) , 

47T 



(48) 



where the function g>(a, 0) is defined by fl4"Dp . By substituting this in ( f26|) and using the fact 
that g(a, (p) is periodic in with period 2n, we obtain 
Now, equation (p9j) can be rewritten as 



/(*) 



1 /" 27r 



32tt 2 
] 

32 



d0 sin%(,2sin#,0)] 2 

I /-27T /"7T /" + 00 

— ■= \ d(p de sine / dipt 5\<p — ip\ + z sin 6* sin ^i) 

^7T 2 JO JO J-00 

/+00 
dip2 S'((p — ip2 + z sin e sin ip 2 ) 
-00 

] /■2'Tr /*7T Z' + OO 

— - / c?0 / de sine / (i'0i 5(0 — ipi + 2 sin e sin-0i) 
Itt 2 Jo Jo j-oo 
r+00 

x / dip2 S"((p — ip2 + z sin e sin -02 ) • 



32tt 



(49) 



The function ip — z sin e sin ip is a monotonously increasing function of ip for |z| < 1, and 
the equation = ip — z sin e sin -0 is satisfied for = by ip = and for = 27r by ip = 2n. 
Hence, this equation for ip has a solution with a value of in [0, 2tt] only if < ip < 2n. 
Hence, after performing the integral, the integration range for ipi becomes [0, 2tt]. Thus, 
we obtain 



/(*) 



1 



2t 



32tt 2 Jo 

where Xi = ipi — z sin e sin ipi, (i = 1,2). Hence, 



de sine / dipt / dip 2 S" (xi — X2) 



(50) 



/(*) 



32vr 
1 

32^ Jo 

v2 



1 r aq ■ f 27T a d ^ [ + °° 
— de sine dxi -t— \ 

i7r z jo jo dxi J-00 



dX2 -1—0 (xi - X2) 
dX2 



de sin e 



2n 



dip 



d 2 ip 
dx 2 



z 



327T 2 



P7T p2tT 

J desin 3 e J dip—— 



dx 
dip 

sin 2 ip 



zsin6 cosip) 1 



(51) 



where we have defined x = ip ~ z sin e sin ip. By integration by parts with respect to ip, we 
obtain 
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where 



h(z) = / d9 sin 2 e / #- — 
Jo Jo (1 — 



COS ^ 



z sin # cos tp) 4 



Now we evaluate this integral. Define 

1 



Hz) 



3 Jo 



d6 sin 6* 



2tt 



di) 



o (1 — ,2 sin 9 cosipy 



The ^ integral can be performed by using the residue theorem after the substitution t 
The result is 



2tt f 71 
H{z) = — d9 sin9 
3 Jo 



1 



3 z 2 sin 2 9 



[l-z 2 sin 2 9f/ 2 2 (1 - z 2 sin 2 9fl 2 



Then, by using the formulae 

dx 



x 



(a + cx 2 ) 3 / 2 a yja + cx 2 



dx 



(a + cx 2 ) 5 / 2 a 2 



x 



cx 



Va + cx 2 3(a + cx 2 ) 3 / 2 



with a = 1 — z 2 and c = z 2 in ( (55|) after the substitution x = cos 9, we have 

4vr 1 



3 (I-2 2 ) 2 

Noting that h(z) = H'(z), we find from this 

16n z 

h{z) = —jr^ 

By substituting this result in (|52| ) we obtain equation (|3 
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FIGURES 



FIGURE 1 




FIG. 1. Scattering potentials Vm and Vs are plotted as functions of r/2M for I = 1, where 
we recall that Vm and Vs are defined as functions of Minkowski and Schwarzschild r coordinates, 
respectively. Asymptotically Vm and Vs fall as 1/r 2 . Vs is only defined outside the black hole 
(r > 2M). Because of the nonexistence of the event horizon in Minkowski spacetime, Vm is also 
defined in the region < r < 2M. 
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FIGURE 2 
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FIG. 2. The radiated power W lr ^ m is plotted as a function of Q for geodesic orbits. Solid and 
dashed lines are associated with the numerical calculation and analytic approximation, respectively. 
The maximum OM considered is 0.068 which is associated with the last stable circular orbit. 
As expected our analytic approximation is only accurate to describe the emission of low-energy 
particles as it can be seen from the coincidence of the numerical and (approximate) analytic curves 
for small fi. It is clear that small angular-momentum waves give the main contribution to the total 
radiated power. 



16 



FIGURE 3 




FIG. 3. The w s > em /W M > em ratio is plotted as a function of as measured by asymptotic 
observers, where the summations involved in the calculation of W s,em and W M,em were performed 
up to I = 3. W s,em and W M,em are the power emitted by an orbiting source as calculated by 
asymptotic observers who assume Schwarzschild spacetime and Minkowski spacetime with Newto- 
nian gravity respectively. The graph is plotted up to Q.M = 0.068, since this is the faster stable 
circular orbit according to General Relativity. Solid and dashed lines are associated with the nu- 
merical calculation and analytic approximation, respectively. Asymptotically (f2 — ► 0) we have 
W S ' em /W M ' em -> 1. As fi increases, however, we see that w S ' em /W M ' em decreases. 
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FIGURE 4 




FIG. 4. Similarly to figure |, the w s > em /W M ' em ratio is plotted as a function of f2. Here, 
however, W M ' em is the power emitted by an orbiting source as calculated by asymptotic observers 
who assume Minkowski spacetime with a gravitational force induced by gravitons, rather than 
with the usual Newtonian force. We note that figure || and figure ^| have similar features: (i) 
w S,emj W M,em _> 1 asymptotically (ft -> 0) and, (ii) w s > em /W M > em decreases up to from 20% 
(see figure |]) to 30% (see figure ||) as Cl increases. 
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FIGURE 5 




FIG. 5. The ratio W s > obs /W S ' em is plotted as a function of f2 for geodesic orbits, where the 
summations involved in the calculation of W s,obs and W S ' em were performed up to / = 3. Solid and 
dashed lines are associated with the numerical calculation and analytic approximation, respectively. 
We note that most of the emitted energy is radiated away to infinity. 



19 



